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We study a d-wave superconductor with dominant d x 2 _ y i -wave order parameter and subdominant pairing 
in either the s- or the d xy -wave channel near a surface. In particular we analyze the influence of surface 
roughness on the mixed order parameter which may break the time-reversal symmetry. We find that the 
subdominant component is suppressed by the roughness independent of its pairing symmetry; for very rough 
surfaces the subdominant component may even vanish completely. Additionally we discuss a possible real- 
valued admixture which counteracts the suppression of the d a .2_„2-wave order parameter at the surface. 



There is an on-going discussion about the surface state of high-T c materials in the superconducting 
phase, which is commonly believed to exhibit d-wave symmetry in the bulk. In particular the question is 
whether a subdominant order parameter occurs at the surface which leads to a breaking of the time-reversal 
symmetry. Theoretical studies show that this is possible for low enough temperatures if the d-wave order 
parameter is strongly suppressed at the surface, and a subdominant pairing interaction exists [l]|2][3|: 
Such a situation is realized, for example, if the surface is oriented in the [110] direction of a d x 2_ y 2- 
wave superconductor with subdominant pairing channels of the s-wave or the d xy -wave type. This would 
manifest itself in the differential conductance of tunnel junctions as a splitting of the zero-bias peak which 
occurs for a tilted single-component d x 2_ y 2 -wave order parameter 1 3 1 . Another experimental possibility is 
to study the magnetic field dependence of the differential conductance, which would also show anomalies 
due to a subdominant pairing interaction |4|. In this work we concentrate on the zero-field case. 

In some experiments on high-T c superconductors a splitting of the zero-bias conductance peak was 
observed 1 5 6 1, which seems to confirm the existence of a subdominant order parameter. On the other hand, 
in numerous similar experiments no such splitting was visible |7]|8l[9). It was shown experimentally for 
YBCO 1 6 1 that one reason for this discrepancy might be a different oxygen doping, i.e. different electronic 
properties of the material. We show here that these results can also be explained by surface roughness, as 
a rough surface can suppress a subdominant admixture of either d^-wave or s-wave symmetry. Moreover 
we demonstrate that in some cases a subdominant admixture to the order parameter occurs without breaking 
the time-reversal symmetry; then the spatial symmetry is modified near the surface. 

To study the superconducting phase of high-T c superconductors we start with a few symmetry con- 
siderations. In conventional superconductors the gauge symmetry is broken only, and usually the order 
parameter is assumed to be isotropic (A(p F ) = A, p^: Fermi momentum). In unconventional super- 
conductors the lattice symmetry is broken in addition to the gauge symmetry, and the order parameter is 
anisotropic. In order to describe high-T c materials we consider a tetragonal lattice symmetry; note that we 
neglect possible orthorhombic distortions, present for example in YBCO. Then the order parameter can be 
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expanded as 

a(pf) = -J2 a Mpf) , (1) 

i 

where rji are basis functions of the irreducible representations of the tetragonal space group, D^h- The 
lowest-order basis functions are given in Tabled Since we neglect the weak dispersion in c-direction, we 
concentrate on three order parameter symmetries, namely d x 2_ y 2-, d xy - and s-wave (see Fig.[Q. In the 
following we will discuss a dominant <i 2 .2_ a 2-wave order parameter, and a weaker pairing interaction in 
either the s-wave or the d x „-wave channel. 

Therefore we consider a superconductor with two (or more) possible order parameter components. Due 
to the pairing interaction each component has its individual critical temperature, T c .i', i.e. in a bulk su- 
perconductor the i th order parameter component would become finite at 7$ if all other components were 
zero. The two cases of interest for us, namely a d x 2_ y 2-wave order parameter with s- or d xy -wave admix- 
ture, are discussed in Refs. 1 10 1 and 1 1 1 1, respectively, with the following results: If both components, the 
dominant and the subdominant, have a finite critical temperature (T c ,i > 0), it is possible that the domi- 
nant order parameter component with T C) i suppresses the subdominant component, which then becomes 
finite at a lower temperature T C 2 < T Cj 2 < T c i, alternatively it may happen that the subdominant pairing 
channel is completely suppressed, i.e. T c a = 0. The latter situation can be assumed for high-X^ supercon- 
ductors, where in the bulk only a d x 2_ y 2-wave order parameter is observed, while other order parameter 
components seem to be absent. 

Although the subdominant pairing can be neglected in the bulk, it might become important near in- 
homogeneities, e.g. impurities or boundaries, which tend to suppress the dominant d x 2_ y 2-wave order 
parameter. At a specular surface, for example, a d x 2_ y 2-wave order parameter which is tilted by 45° 
with respect to the surface normal is suppressed completely. In the following we will discuss the inter- 
play between the suppression of the d x 2_ y 2-wa\e order parameter near an interface and the existence of a 
subdominant order parameter component in this region. In particular we study the influence of interface 
roughness, which also modifies the behavior of the d x 2_ y 2-wa\e component. 



Dih'- basis functions 

Ai g Vi{p) = + PyiPi (anisotropic) s-wave 

Mg T]2 (p) = PxPy (pi ~ P% ) 

Big V3(p) =P 2 X -P 2 y d x 2_ y 2-wave 

B 2g V4(p) =PxP v d xy -wave 

E g %(P) = (PxPz,PyPz) 

Table 1 Basis functions of lowest order for the group D^. 




Fig. 1 Angular dependence of the order parameter for three different symmetries: s-, d x 2_ y 2-, and d xy - 
wave. 
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Here we apply the theory of quasiclassical Green's functions in thermal equilibrium 1121 1131 . The 
Green's functions are determined by the Eilenberger equation, 

T 3 E + eT 3 v F -A(r)+iA(p F ,r),g(E,p F ;r) + zv F • d r g(E, p F ; r) = , (2) 

supplemented by the normalization condition 

[g(E,p F ;r)] 2 = i. (3) 

Here ?i are the Pauli-matrices in Nambu-space. The (spin-singlet) order parameter 

( A .£, r) ^r > ) 

must be determined self-consistently via 

A(p F ,r) = -ttA/'oT ]T {V(p F ,p' F )g(iE n ,p' F ;r)) p , p , (5) 

\E n \<E a 

where (. . .} , denotes an average over the Fermi surface (assumed to be spherical). The cut-off energy is 
E c , and the normal-state density of states at the Fermi energy is denoted by Mo per spin. According to the 
above remarks on the order parameter symmetry we can expand the pairing interaction as follows: 

V(p F ,p' F ) = -^2ViT)i[p F )r]i(p F ) . (6) 

i 

The current density can be determined from the Greeen's function via 
j(r) = -ienMoTj2^ [** (vF^K,PF;r)) pp 

E, 



(7) 



and the angle-resolved density of states is given by 

M(E, p F : r) = -y Re{Tr [f 3 g(E + t0 + , p F ; r)] }. (8) 
The density of states at the surface is connected to the differential tunnel conductance by the relation 1141 

G(V) = Ae 2 (v F . x T(p F )M(eV, p F ; 0)) pp , (9) 

where T(p F ) = Top Fx /p F (% -C 1) is the angle-dependent transparency, and A is the cross-section of 
the tunneling contact; this results in the normal state resistance i?^ 1 = 4:Ae 2 Mov F To/3n. 

The quasiclassical approach is valid only on scales which are large compared to the Fermi wavelength, 
and is not directly applicable near interfaces. This problem must be circumvented using boundary condi- 
tions for the Green's functions. For specular interfaces, i.e. when the momentum parallel to the interface 
is conserved, one can use Zaitsev's boundary conditions 1 15 1; these can be extended to include interface 
roughness 116lll7lfT8lll9l . We use here an approach which was introduced by Shelankov and Ozana |20|, 
in which the properties of the interface are described by a scattering matrix. 

In order to describe interface roughness, which is present on the atomic scale, we choose the scattering 
matrix to be a random unitary matrix 1161 1141 12T1 . Here the only relevant roughness parameter is the 
probability for specular scattering, \u\ 2 . A rough interface leads to a reduced weight for specular scattering, 
\u\ 2 < 1, and with a probability of 1 — \u\ 2 a particle is scattered randomly into any other direction. This 
model is described in detail in 1141 1211 . The disorder averaged quantities are denoted by (. . . ). 



4 



T. Luck: d-wave superconductors near boundaries 




The superconductor is described by a dominant d x 2_ y 2 -wave order parameter with a critical temperature 
T c , and a second order parameter component of either s- or d^-wave type with an individual critical 
temperature T c .2 = 0.3T C ; this choice ensures that in the bulk only a d 2 .2_ y 2-wave order parameter is 
present. The surface normal is in the x-direction, and we assume translational invariance in the y-direction 
(parallel to the surface). 

The behavior of a single-component d x 2_ y 2-wave superconductor near a surface is well known (1611171 
1181 1191 fl4l . The following results are of particular importance for our purpose: If the order parameter 
is tilted with respect to the surface normal, trajectories exist where the order parameter changes its sign 
due to surface scattering; this leads to a zero-bias peak in the differential conductance of tunnel contacts, 
and a suppression of the order parameter near the surface. Without surface roughness and for a tilting 
angle of 45° the order parameter at the surface is zero. Surface disorder leads to a finite value of the order 
parameter, and in most models to a broadened zero-bias conductance peak. Furthermore, even an untilted 
d x 2_ y 2-wave order parameter is suppressed near a rough surface 1141 . 

First we focus on a (4,2 _ y 2 -wave superconductor with a subdominant pairing in the s-wave channel; the 
order parameter is rotated by 45° with respect to the surface normal. Here it is important to note that an 
isotropic s-wave order parameter itself is inert against scattering at a rough surface. Without roughness 
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(\u\ 2 — 1) the dominating d x 2_ y 2-v/a\e component is completely suppressed at the surface (x = 0), and 
a finite s-wave admixture occurs near the boundary with a relative phase ip = ±7r/2 (Fig. O. This state 
has a spontaneously broken time-reversal symmetry near the surface. As a consequence a current flowing 
in y-direction is present (Fig. |2J, and a splitting of the zero-bias conductance peak is found (Fig.0}. These 
are well-known results for the clean case p|l2l|3l l22l . 

For a rough surface the ^2 -wave component recovers and has a finite value at the surface; this 
in turn leads to a suppression of the subdominant component, the suppression becoming stronger with 
increasing roughness (Fig.|2ji. As a consequence the splitting of the zero-bias conductance peak and the 
parallel current are reduced. For very strong roughness (\u\ 2 < 0.08) no subdominant order parameter is 
observed. Therefore the differential conductivity is no longer affected by the subdominant pairing, and no 
parallel current exists. 

Next we consider the situation with a subdominant d^-wave component. In the clean case we find 
similar results as above, since the 45°-rotated d xy -wa\e order parameter is not suppressed by the surface. 
Therefore a finite d xy -wa\e admixture with a relative phase of ip = ±7r/2 exists near the surface, and the 
resulting state exhibits broken time-reversal symmetry. Consequently a current in y-direction flows near 
the surface (Fig.|3j, and the zero-bias conductance peak splits (Fig.|3. 
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Here the current behaves qualitatively different (Fig. to what we found for subdominant s-wave 
pairing (Fig. [3]»; in particular here the current changes its sign at some distance from the surface. This is 
related to the fact that a mixed d x 2_ y 2 + id xy order parameter has a finite magnetic moment, which is not 
the case for a d x 2_ y 2 + is order parameter. 

A crucial difference occurs for rough surfaces: The d^-wave component is not only suppressed by the 
dominant order parameter but, as it is anisotropic, also by the surface roughness. Therefore the subdomi- 
nant admixture decreases rapidly with increasing surface roughness; we find a full suppression already for 
M 2 < 0.42 (Fig.©. 

We now study the situation when the order parameter is rotated by only 24° with respect to the surface 
normal, and a subdominant d xy -wa\e pairing is present. The crucial difference to the previous cases is that 
the dominant order parameter is not completely suppressed by the specular surface. We also find a finite 
value of the subdominant order parameter but with a relative phase tp « 0; i.e. the admixture is mainly 
real-valued, and the time-reversal symmetry-breaking imaginary part is negligible (Fig. Therefore 
no splitting of the zero-bias conductance peak can be observed (Fig. |6j. The real-valued subdominant 
component leads to an effective rotation of the dominant order parameter such that the tilting angle is 
diminished. As a result the suppression of the order parameter at the surface is reduced; the differential 
conductance shows that spectral weight is shifted from lower to higher energies compared to the case 
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Fig. 5 Dominant d a .2_ !/ 2-wave order parameter with (real) subdominant d xy -wave component, for a = 
24° atT = 0.1T C ; A = 2.14T C , and Co = v f /ttA . 
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Fig. 6 Differential conductance versus voltage for a = 24°, for subdominant pairing in the d^-wave 
channel. Note, that for the clean case (|u| 2 = 1) the zero-bias conductance peak has been broadend by 
adding a finite imaginary part to the energy. 



without subdominant admixture. For subdominant s-wave pairing a similar effect occurs. For a 45°-tilting 
of the order parameter no real-valued admixture develops due to the symmetry of the system. 

In conclusion we have shown that, for a 45°-rotated order parameter, surface roughness suppresses 
the subdominant admixture in the d^-wave as well as in the s-wave case. We have also found that the 
subdominant component may vanish completely for very strong roughness. As expected the cZ^y-wave 
component is more strongly affected by roughness than an s-wave component. For a tilting angle a < 
45°, a real-valued admixture occurs as well; the order parameter is modified so that its suppression is 
weakened. Our results imply that surface roughness plays a crucial role for the experimental observation 
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of subdominant pairing. In particular the roughness must be small enough to detect the splitting of the 
zero-bias conductance peak. 
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